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ABSTRACT 

We have developed a method for measuring higher-order weak lensing distor- 
tions of faint background galaxies, namely the weak gravitational flexion, by fully 
extending the Kaiser, Squires & Broadhurst method to include higher-order lens- 
ing image characteristics (HOLICs) introduced by Okura, Umetsu, & Futamase. 
We take into account explicitly the weight function in calculations of noisy shape 
moments and the effect of higher-order PSF anisotropy, as well as isotropic PSF 
smearing. Our HOLICs formalism allows accurate measurements of flexion from 
practical observational data in the presence of non-circular, anisotropic PSF. 
We test our method using mock observations of simulated galaxy images and 
actual, ground-based Subaru observations of the massive galaxy cluster A1689 
{z = 0.183). From the high-precision measurements of spin-1 first flexion, we 
obtain a high-resolution mass map in the central region of A1689. The recon- 
structed mass map shows a bimodal feature in the central 4' x 4' region of the 
cluster. The major, pronounced peak is associated with the brightest cluster 
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galaxy and central cluster members, while the secondary mass peak is associated 
with a local concentration of bright galaxies. The refined, high-resolution mass 
map of A1689 demonstrates the power of the generalized weak lensing analysis 
techniques for quantitative and accurate measurements of the weak gravitational 
lensing signal. 

Subject headings: cosmology: theory — dark matter — galaxies: clusters: indi- 
vidual (A1689) — gravitational lensing 

1. Introduction 

Propagation of light rays from a distant source to the observer is governed by the 
gravitational field of intervening mass fiuctuations as well as by the global geometry of the 
universe. The images of background sources hence carry the imprint of the gravitational 
potential of intervening cosmic structures, and their statistical properties can be used to test 
the background cosmological models. 

Weak gravitational lensing is responsible for the weak shape-distortion and magnification 
of the images of background sources due to the gravitational field of intervening matter (e.g., 
Bartelmann & Schneider 2001; Umetsu, Futamase, & Tada 1999). To the first order, weak 
lensing gives rise to a few - 10% levels of elliptical distortions in images of background sources, 
responsible for the second-order derivatives of the gravitational lensing potential. Thus, 
the weak lensing signal, measured from tiny but coherent quadrupole distortions in galaxy 
shapes, can provide a direct measure of the projected mass distribution of cosmic structures. 
However, practical weak lensing observations subject to the effects of atmospheric seeing, 
isotropic/anisotropic PSF, and (residual) camera distortion across the field of view, which 
must be examined from the stellar shape measurements and corrected for in the weak lensing 
analysis. Practical methods for PSF corrections and shear measurements/calibrations have 
been studied and developed by many authors, such as pioneering work of Kaiser, Squires, & 
Broadhurst (1995, hereafter KSB), the Shapelets technique which describes PSF and object 
images in terms of Gaussian-Hermite expansions (Refregier 2003), and recent systematic, 
collaborative efforts by The Shear TEsting Programme (Heymans et al. 2006; Massey et al. 
2007). 

Thanks to these successful developments in weak lensing techniques as well as in instru- 
ment technology, the quadrupole weak lensing has become one of the most important tools 
in observational cosmology to map the mass distribution in individual clusters of galaxies 
(e.g.. Kaiser & Squires 1993; Broadhurst et al. 2005a; Okabe & Umetsu 2007; Umetsu & 
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Broadhurst 2007), measure ensemble-averaged mass profiles of galaxy-group sized halos from 
the galaxy-galaxy lensing signal (e.g., Hoekstra et al. 2001; Hoekstra et al. 2004; Parker et 
al. 2005; Mandelbaum et al. 2006), study the statistical properties of the large scale struc- 
ture of the universe from the cosmic shear statistics (e.g.. Bacon et al. 2000; van Waerbeke 
et al. 2001; Hamana et al. 2003), and search for galaxy clusters by their mass properties 
(e.g., Schneider 1996; Erben et al. 2000; Umetsu & Futamase 2000; Wittman et al. 2001). 

In recent years, there have been theoretical efforts to include the next higher order 

distortion effects as well as the usual quadrupole distortion effect in the weak lensing analysis 
(Goldberg & Natarajan 2002; Goldberg & Bacon 2005; Bacon et al. 2006; Irwin & Shmakova 
2006; Goldberg & Leonard 2007; Okura, Umetsu, & Futamase 2007). We have proposed in 
Okura, Umetsu, & Futamase (2007, hereafter OUF) to use certain convenient combinations 
of octopole/higher multipole moments of background images which we call the Higher Order 
Lensing Image's Characteristics (HOLICs), and have shown that HOLICs serve as a direct 
measure for the next higher-order weak lensing effect, or the gravitational flexion (Goldberg 
& Bacon 2005) and that the use of HOLICs in addition to the quadrupole shape distortions 
can improve the accuracy and resolution of weak lensing mass reconstructions based on 
simulated observations. 

Recently, Goldberg & Leonard (2007) extended the HOLICs approach for flexion mea- 
surements to include observational effects, namely the Gaussian weighting in shape-moment 

calculations (see the appendix therein) and the isotropic PSF effect, under the assumption 
that PSF is nearly circular, and tested their extended HOLICs approach with simulated and 
HST/ACS observations. Leonard et al. (2007) have applied the extended HOLICs method 
to reconstruct the projected mass distribution in the central region of the massive galaxy 
cluster A1689 a.t z = 0.183, and revealed substructures associated with small clumps of galax- 
ies. Further Leonard et al. (2007) found that in dense systems such as galaxy clusters the 
HOLICs technique is robust and less sensitive than the Shapelet technique to contamination 
by light from the extended wings of lens/foreground galaxies. 

In the present paper we develop a method for measuring flexion by the HOLICs approach 
by fully extending the KSB formalism; We take into account explicitly the effects of Gaussian 
weighting in calculation of noisy shape moments and higher-order PSF anisotropy as well 
as isotropic PSF smearing. We then apply our method to actual, ground-based Subaru 
observations of A1689, and perform a mass reconstruction in the central region of A1689. 

The paper is organized as follows. We first summarize in §2 the basis of weak gravita- 
tional lensing and the flexion formalism. In §3, we derive the relationship between HOLICs 
and flexion by incorporating Gaussian smoothing in shape measurements in the presence of 
isotropic and anisotropic PSF. The practical method to correct the isotropic/anisotropic PSF 
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effects will be presented in §4. In Section 5 we use simulations to test our flexion analysis 
method based on the HOLICs moment approach. We then perform a weak lensing flexion 
analysis of A1689 by our fully-extended HOLICs approach, and perform a mass reconstruc- 
tion of A1689 from the HOLICs estimates of flexion. Finally summary and discussions are 
given in §6. We refer interested readers to a complete appendix^ for details of the derivation 
of flexion-observable relationships in practical observations. 



2. Basis of Weak Lensing and Flexion 

In this section, we summarize general aspects of weak gravitational lensing and flexion 
formalism, following the complex derivative notation developed by Bacon et al. (2006). A 
general review of quadrupole weak lensing can be found in Bartelmann & Schneider (2001). 



2.1. Spin Properties 

We define the spin for weak-lensing quantities in the following way: A quantity is said 
to have spin if it has the same value after rotation by 2tt/N . The product of spin- A and 
spin-i? quantities has spin {A + B), and the product of spin-A and spin-_B* quantities has 
spin {A — B). Then, as we shall see in the next subsection, the lensing convergence k is 
a spin-0 (scalar) quantity. The complex shear 7 and the reduced shear g = 7/(1 — k) are 
spin-2 quantities. The first and second flexion fields, F and G, are a spin-1 and a spin-3 
quantity, respectively. 



2.2. Weak Lensing and Flexion Formalism 

The gravitational deflection of light rays can be described by the lens equation, 

/3 = e-V^(0), (1) 

where il^{6) is the effective lensing potential, which is defined by the two-dimensional Poisson 
equation as V'^ip{6) = 2k{6), with the lensing convergence. Here the convergence k = 
J dTim^cvlt is dimensionless surface mass density projected on the sky, normalized with 
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respect to the critical surface mass density of gravitational lensing, 

Ds 

where D^, Dg, and are the angular diameter distances from the observer to the deflector, 
from the observer to the source, and from the deflector to the source, respectively. By 
introducing the complex gradient operator, d = 81+182 that transforms as a vector, 8' = 86"^^, 
with <f> being the angle of rotation, the lensing convergence k is expressed as 

K = ^88*ij, (3) 

where * denotes the complex conjugate. Similarly, the complex gravitational shear of spin-2 
is defined as ^ 

7 = 7i + Z72 = -88iIj. (4) 

The third-order derivatives of iplO) can be combined to from a pair of the complex flexion 
fields as (Bacon et al. 2006): 

= T^ + iT2 = (5) 

Q = Qi+iQ2 = \888ij. (6) 

If the angular size of an image is small compared to the scale over which the lens potential 
il) varies, then we can locally expand the lens equation ([T]) to have: 

d(3i = Aijddj + ^Vijkd9jd9k (7) 

to the second order, where Aij is the Jacobian matrix of the lens equation and Vijk = Aij^k = 
—ip^ijk is the third-order lensing tensor. 



A- 



1 - K - 7i -72 
-72 1 - K + 71 

T-^ijk ■^ijk ~\~ Qijk- (9) 



The third-order tensor 'Dijk can be expressed with the sum of the two terms, Pjjfc = J-'ijk+Qijk, 
with the spin-1 part J-'ijk and the spin-3 part Qtjk, composed of the real/imaginary part of 
the flexion flelds: 
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Note that flexion has a dimension of inverse length (or inverse angle), meaning that the 
flexion effect depends on the angular size of the source. The shape quantities affected by 
the first flexion T alone have spin-1 properties, while those affected by the second flexion 
Q alone have spin-3 properties. These third-order lensing fields naturally appear in the 
transformation equations of HOLICs between the lens and source planes. 



2.3. Quadrupole Lensing Observable - Complex Ellipticity 

In the KSB approach, we use quadrupole moments Qij of the surface brightness distri- 
bution of background images for quantifying the shape of the images: 

JdWqj[Iie)] ' ^ ' 

where qi[I{0)] denotes the weight function used in noisy shape measurements and A9i = 
6i — 6i is the offset vector from the image centroid. The complex ellipticity x is then defined 
as 

^ ^ Qu - Q22 + 2iQi2 ^^2^1 

Qii + Q22 

The X transforms under the lens mapping as 

^ l + \g\^-2Re[gx*Y ^ ^ 

where (7 = 7/(1 — k) is the spin-2 reduced shear. In the weak lensing limit (k, I7I ^ 1), 
equation (IT^ reduces to x^^^ ~ X ~ 27. Assuming the random orientation of the background 
sources, we average observed ellipticities over a sufficient number of images to obtain 

(X) -2g^ 27. (15) 



2.4. Centroid Shift due to Lensing 

In the moment methods such as KSB, we quantify the shape of an image by measuring 
various moments of I {6), in which the moments are calculated with respect to the observable 
centroid of the image 6, or the center of the light, defined by the first moment of I{0). 
However, in the presence of gravitational lensing, this apparent center can be different from 
the point 0{p) that is mapped using the lens equation ([T]) from the center of the unlensed 
light, ^. We refer to this point as the "true" center of the image. The difference between 
these two centers, namely, the apparent and the true centers, causes a significant effect in 
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evaluating the first flexion, as pointed out by Goldberg & Bacon (2005). The relationship 
between the two centers is 

em « e, - trQ Qf, + ^[F*x]i + l[Gx*l^ =0i- Ai,„ (16) 

where trQ = Qu + Q22 is the trace of Qij, F and G are the reduce Flexion, defined by 
F — J-'/ {I — k) and G = Q/{1 — n), respectively, and A^^i is the displacement vector from 
the true to the apparent center due to gravitational lensing. Therefore, by taking into account 
the centroid shift in shape measurements, we have the relation between A/3j — (5i — and 
A6li ^Oi-Oi as 

^ AijAOj + -VijkAejAOk + trQ ( -F, + -[F*x]i + -^[Gx*]i] 

= AijAej + ^VijkAejAek + AL,i. (i7) 

A detailed derivation of the above relation is given in Appendix B of OUF. 



2.5. Flexion Observable - HOLICs 

The ficxion fields, JF and Q, can be measured from proper combinations of higher-order 
shape moments with the corresponding spin properties and the dimension, as explicitly 
shown by OUF. Higher-order moments of images are defined as a straightforward extension 
of the quadrupole moment. The octopole moment Qijk and the 16pole moment Qijki are 
defined as follows: 

^ jd'9qj[i{d)]AeAeA0k .... 
^''^ - jd?eqi[i{e)] ^ ' 

^ jd?eqi[i{e)\AeiAejAekAei 
- jd?eqi[i{d)\ ■ ^ ' 

Then, C, and 5 of the spin-1 and spin-3 HOLICs, respectively, are defined by 

Qui + Q122 + i (Q112 + Q222) 



Qui ~ 3Q122 + i (3Q112 — Q222) ^2Q"^ 



where ^ is the spin-0 normalization factor. 



C — Qllll + 2Q1122 + Q2222- 



(21) 



Finally, the transformation equations between unlensed and lensed HOLICs are obtained 
as (see OUF) 

C - 2^C* - 9*^ - ii8F*r] - lQ^^^F*x + 9F- U^^^F + 2Gr]* - 2^^Gx* + G*X) 
^ ~ (1 - k)(1 - 4Re [g*r]] - 5Re [Fl*] - Re [Gi}jj]) ' 

6-3gC-jU0Fr] + 7F*X-18^^Fx + SG) 

{1- K){l-4Re [g*r]]-5Re[FL*j] -Re [Gl}jj])' ^ ' 

where t] and A are dimensionless spin-2 and spin-4 quantities, respectively, defined with 
16-pole moments, and t/, Lm and Ly are spin-1, spin-3, and spin-5 quantities, respectively, 
defined with 32-pole moments (see Appendix and OUF). We note that, the above equations 
( I22l) and (!22l) are obtained under the sub-critical lensing condition, i.e., det^ > (see 
Schneider & Er 2007). 

Since the HOLICs C and 5 are non-zero spin quantities with a direction dependence, 
the expectation value of the intrinsic C, and 5 are assumed to vanish. To the first order in 
fiexion, we have the linear relations between the HOLICs and fiexion fields as 

^ ^ (23) 



9 _ Q [trQf 
4 '3 5 



G ^ ^{6). (24) 



3. HOLICs and Flexion in Practical Applications 

For a practical application of the HOLICs approach, we must take into account various 
observational effects such as noise in the shape measurement due to readout and/or sky 
background and the dilution of the lensing signal due to the isotropic/anisotropic PSF effects. 
Thus, one cannot simply use equations fl22l) and fl22|) to measure the flexion flelds. In this 
section, we introduce the Gaussian weighting in moment calculations, as done in the KSB 
formalism for quadrupole weak lensing, and derive the relevant transformation equations 
between unlensed and lensed HOLICs by taking into account explicitly the effect of Gaussian 
smoothing. 



3.1. Redefining HOLICs for Noisy Observations 



Now we introduce a weight function lV(|A^p/(T^) having a characteristic width a for 
practical, noisy moment measurements, redeflning the octopole and 16-pole moments of the 
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brightness distribution I{d) as 

Q^Jk = J d^eiie)AeiAejAekW{\Ae\ya^), (25) 
Q^jkl = j Sei{e)Ae,AejAekA9iW{\Ae\^ /a^), (26) 

where Qs here are no longer normahzed with the corresponding flux of the image (see equation 
[12]). The redefined shape moments enter equation (1201) . We provide in Appendix detailed 
definitions of HOLICs in practical applications. 

3.2. Lensing-Induced Centroid Shift in Weighted Moment Calculations 

When weighted moments are used for calculating HOLICs, the centroid shift due to 
lensing ffTUj) . to the first order, is changed in the following way: 

3trQ° _|_ 3g°' 



AL = e- eiP) ^ ^ A^F, (27) 



where 9 = 9i + 162 and 6{J3) = 6i{J3) + 26*2 (/?) are the apparent and the true centers of 
the image (see §2.4p . respectively, in the complex form calculated using the weight function 
W{x), M = J (P 9 1 {9)W (lAOl"^ / a^) is the monopole shape moment (or flux), quantities with 
subscript "a" represent those calculated with respect to the apparent center, and quantities 
with prime represent those measured with W'{x) = dW{x) / dx as the weight function instead 
of W{x)] A° is the spin-0 coefficient in Al. The deviation from unity in the denominator of 
A° is obtained by properly expanding the weight function W(x) in moment calculations, and 
this term does not appear in Goldberg & Leonard (2006) 's formulation. Hence, the complex 
displacement from the true image center, A6^, can be expressed in terms of that from the 
apparent image center, A^*^, and the complex centroid shift, Al, as 

A^* = Ar + Al. (28) 



For interested readers, we refer to Full Appendix B.l for detailed calculations of the lensing- 
induced centroid shift with a weight function. 
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3.3. Relation between the Weighted HOLICs and Flexion 

In weighted moment calculations, the transformation equations between HOLICs and 
flexion must be modifled accordingly. To the flrst order, we have 

^0 = Qllllll + 3(5ii2222 + 3(5iiii22 + Q2222225 (31) 

where quantities with subscript "t" represent those calculated with respect to the true center, 
or using A6'*. By using equation fpH|) . we can express and 5* in terms of practically 
observable quantities (with subscript " 

- + 2^A. + i^A. = C+ + i^Ao) F, (32) 

5* ^ 5" (33) 

to the flrst order. Here the term -^^Al in equation ( l32l) is again caused by the centroid 
shift in the weight function (a similar, but different, expression was obtained by Goldberg 
& Leonard 2007). Finally, we obtain the following transformation equations in the case of 
weighted moment calculations: 



Sis) 



1- K 
1 



f9 3< trQ'^ e' .oV' 



1-K 

We show the details of these calculations in Full Appendix B.3. 



(34) 
(35) 



4. Isotropic and Anisotropic PSF Corrections for HOLICs Measurements 

In this section we present a detailed prescription for the PSF anisotropy and circular- 
ization correction in HOLICs-based flexion measurements by extending the KSB formalism. 
We closely follow the treatment and the notation given in §4.6.1 of Bartelmann & Schneider 
(2001). 
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4.1. General Description for PSF 

The observed surface brightness distribution I°^^{0) can be expressed as the true surface 
brightness I{d) convolved with an effective PSF P{6), 

r^'{e)^ j (f^i{^)p{o --a). (36) 

Following the KSB formalism, we assume that P is nearly isotropic, so that the anisotropic 
part of P is small. We then define the isotropic part of P as the azimuthal average over P, 
and decompose P into an isotropic part, P*^", and an anisotropic part, g, as 

P(t?) = J(f<l>q(cl>)P'^\^-cf>), (37) 

where both P**° and q are normalized to unity. We then define P^°{d) as the surface bright- 
ness distribution smeared by the isotropic part P**", 

f'°{e)^ j (f^i{e)p''°{e-'&). (38) 

The observed surface brightness P^'^id) is obtained by convolving P'^°[0) with the anisotropy 
kernel q{d) as 



Tobs 



(0) = J (P^q{e - '0)P'"{'0). (39) 



In the original KSB method for quadrupole weak lensing, only the spin-2 PSF anisotropy 
described by the quadrupole moment of the anisotropy kernel q has to be taken into account. 
However, in order to correct HOLlCs with spin-1 and spin-3 properties for the anisotropic 
PSF effects, we need to take into account the corresponding dipole and octopole moments 
of the anisotropy kernel q having spin-1 and spin-3 properties. We expand the integral of 
arbitrary function f{0) and I°^^{0) to obtain 

£ef{d)F'%d)^ j£cj>p^°{4>) jd'9f{d + 4>)q{e)^ jd'cj>p^°{4>)f{4>) 

+q, jd'<t>P'\4>)^^+\qki Jd'^P^%4>)- 



d(f)kd(f)i 



+ ^qkim I d'</^r%cP) + , (40) 

6 J d(j)kd(j)id(j)m 



where 



q,^ jd^<t>q{(t>)<t>i, (41) 
Qij ^ I(P<p(l{(l>)<Pi<f>j, (42) 
qijk = Jd^(j)q{(l))(j)i(j)j<f)k (43) 
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are the dipole, quadrupole, and octopole moments of the PSF anisotropy kernel q, respec- 
tively. 



4.2. Centroid Shift due to PSF Anisotropy 

The PSF anisotropy can cause a centroid shift between the images defined in terms of 
P^°{d) and I°^^{d). This centroid shift is essentially due to the spin-1 PSF anisotropy. 

The observed weighted-center of image can be expressed in the complex form as 

jdwi°bs(e)wi\e\ya^) ' ^ ' 

where 6 = 6i + 162 is the complex angular position. Similarly, the weighted-center of a 
hypothetical image defined in terms of /**° is 

j<pei''°{6)w{\e\'/a^) ■ ^ ' 

Then, the offset between the two centers 0°^^ and 9'^^° is given as 

Ap = 9°^' - 9''° = A^^^" - M"^', (46) 

where A^'*° = 9 — 9^^° and A^"*"* = 9 — 9°^^ are the complex displacements to an arbitrary 
point 9. Expanding Ap to the first order of q yields 

Ap = A^-° - A9'^'^ ^D, + M!^^_^[^t_^^Q ^ d,+ -f C„ (47) 

where quantities with // and /// refer to those calculated with W"{x) and W"'{x) as the weight 
function, respectively, {Dg, (g) are spin-1 complex moments of the anisotropy kernel defined 

by 



D, = J d'Him (48) 
C, = f d'M<P)<P<P<P*, (49) 



(p = 0i + i02, (50) 

and {P^, P^) are certain combinations of shape moments defined as the coefficients in equa- 
tion (Cll) of Full Appendix C.1.2, associated with the complex dipole moment D with spin-1 
properties. Wc note that Ap, the PSF-induced centroid shift, cannot be constrained from 
observations; however, as we shall see in the next subsection, one can constrain the spin-1 
octopole component in Ap, namely (g oc Ap — Dg, in the first order approximation. More 
detailed calculations are presented in Full Appendix C.1.2. 
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4.3. Anisotropic PSF Correction for HOLICs 

We use equation (jlU]) to relate observable HOLICs (with subscript "obs") to those 
defined in terms of /*^° (with subscript "iso"). To the first order of g, we have the following 
relations for the spin-1 and spin-3 quantities: 

r 



D 



C'' ^ r° + U2trQ + ^VD,-Ap) + UM+^+^^" ■ 



C° + (P'c - ^Pt) C„ (52) 



5"'' - ^^^° + ^(^ + ^ + S + Sj^^ = ^^^° + ^°^^' (53) 
where we have deffined the weighted first moment D by 

D = ^ ^'^'^ ' , (54) 

and the spin-3 PSF anisotropy 6q by 

Sg= [ d^e q{(l))^(l)^. (55) 



In general, the PSF varies spatially over the field. If the spatial variation of PSF is sufficiently 
smooth, then one can measure (q and 6q for a set of stars. Since C^^" and 6^^° vanish for stars, 
the spin-1 and spin-3 PSF anisotropies can be obtained as 

C. = . / .0^* .X , (56) 



f pO _ Pd pA 



= (57) 

where quantities with asterisk denote those measured for stellar objects. Note that unlike 
the higher-order PSF anisotropies (q and 6q, Ap cannot be determined from observations, so 
that -0°^'' cannot be corrected for the anisotropic PSF effect. We show detailed calculations 
in Full Appendix C.1.3. 



4.4. Isotropic PSF Correction for HOLICs 



This subsection provides the method of correcting for the isotropic PSF effect on the 
flexion measurement. Firstly, from Liouvelle's theorem we have 1(6) = I^^\f3{6)) with 
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being the surface brightness distribution of the unlensed source. We then consider 

= [ d'^§&\^)P''°i0 - A-^^) = i(Ae), (58) 

det^ J 

where we have defined the brightness distribution convolved with the hypothetical PSF 

i{/3) = jd^I^'\<t>)P{(3-<t>). (59) 

Pm = (60) 

The P can be regarded as an effective PSF relating / to I^'^^; the anisotropic part of P is 
caused by lensing. The octopole moment defined with /(/?) is written as 

Qijk = J d''(3Al3,Al3,ApJ{f3)W • (61) 

Then, ( and 6 of HOLICs are defined in terms of Qijk and Qijki as 

P Qui + Ql22 + iiQlU + Q222) 

C = , (62) 

Qiiii + 2Q1122 + Q2222 

^ _ <5iii — 3(5i22 + ii^Qiu — Q222) j-gg.^ 

Qllll + 2(5ll22 + Q2222 

Substituting the expression for A/5 into the above equations and using equation ( IMl) . we 
obtain 

\ + ^- - ^,Kj F = C"° - (C» + Cf A» ) F, (64) 

S « r° - + G = CjG (65) 

to the first order. Here C°s are certain combinations of shape moments defined in terms 
of (see Full Appendix C.2). In practice, however, one can replace J*'^° with 1°^^ for 
calculating C°s of first order in flexion. 

Next, we decompose P{/3) into an isotropic part, P*'^°, and an anisotropic part, q, as 

P{(3) = jd^m-4>)P''\4>)- (66) 
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With P*^° we define the surface brightness distribution j° smeared by the isotropic part P* 



as 



(67) 



Then, the relationship between /(/3) and /°(/3) is the same as that between I"'''^{0) and 
P'^°{6) except that P instead of P. Therefore, we obtain the following relations to the first 
order: 



c 

6 



5' + PK 



where G and 6q are defined as 



£9 



(68) 
(69) 

(70) 
(71) 



P^ and P^ are spin-0 coefficients defined by certain combinations of shape moments (see Full 
Appendix C.2), and are calculated using J° instead of P'^" to the first order approximation. 
For stellar objects, C° = 5° = 0, so that higher-order PSF anisotropics of q are obtained 
using equations ( l64l) and ( l68l) as 



P9 - &P^ 



Further, to the first order approximation, we have 



-G. 



P^ 



(72) 
(73) 

(74) 
(75) 



Finally, we can relate unlensed HOLlCs in terms of P{/3) to observed HOLlCs corrected for 
the PSF anisotropy by 



5° 



1 - n 



(C^j-C^A^) / po 
C - {Cc + A,) P + ^—-j,^-^ \^P^ - j^P^ \ F 



(f<0\ 

Vs )* 



(76) 
(77) 
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Assuming {(^) = ((5°) = for unlensed sources, we obtain the desired expressions for flexion 
as 




A detailed derivation of the above equations is provided in Full Appendix C.2. 



5. Simulations and Observations 

5.1. Simulated PSF Anisotropies and Corrections 

We use simulations to test and assess the limitations of our PSF correction scheme 
for the flexion measurement. To do this, we assume particular models for describing the 
isotropic/ anisotropic PSF and the surface brightness distribution for a source. In the present 
simulations observational noise and lensing effects are not taken into account. 

First, we assume for the stellar surface brightness distribution a two-dimensional Dirac 
delta function, h{0) — Sjj{0), while for the galaxy surface brightness distribution a truncated 
Gaussian as 



1 


\e\ 




1 - exp 1 


/ d2 \ 
-^max j 







1,^,(0) = exp - exp --^ for \e\ < R^^, (80) 



where i?gai and -Rmax are the Gaussian dispersion and the truncation radius of /gai, respec- 
tively. In the following we set /?max = 3i?gai. Next, we assume the isotropic part of PSF, 
P*^°, also follows a truncated Gaussian of the form: 



(0) = ^ exp -i-^ - -— ^ exp --^ ) for \6\ < 9^^, (81) 



where aiso and ^max are the Gaussian dispersion and the truncation radius of P**°, respec- 
tively. In the following we set O^ax — ScTjso. Finally, we adopt the PSF anisotropy kernel 
q{6) of the following form: 

q{d) = Aaniso 7^ for |6>| < ^aniso, (82) 
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where Aaniso is the normalization factor that controls the strength of PSF anisotropy, and 
^aniso is the truncation radius. We have chosen the direction of anisotropy along the x-axis 
(6'i-axis). Thus, the anisotropy kernel q{6) has two free parameters, (Aaniso, 6'aniso)- We vary 
the parameters (Aaniso, 6'aniso) to test our anisotropic PSF correction scheme as a function of 
degree of PSF anisotropy. We take 6'aniso in the range of 6'aniso ^ (0, (Jiso)- 

Having set up the models, we then produce pixelized images for model stars and galaxies 
using the surface brightness distributions h{6) and /gai(^), which are then convolved with 
the model PSF to yield If '{6) and No observational noise or intrinsic/gravitational 

flexion has been added in the present simulations. Here we consider the following set of 
Gaussian source radii, -Rgai = IcTj^o, 2crjso, 3o"jso- For each set of (Aaniso, 6'aniso) , we measure 
various PSF moments, such as Dg — Ap oc (q, from mock stellar images using the Gaussian 
weight function of dispersion Vg = aiso- On the other hand, we measure various shape 



moments for mock galaxy images of I^\{d) using the Gaussian weight of = y-Rgai + (^ho 
(i.e., Gaussian dispersion of the PSF-convolved image). Then, we correct observed HOLICs 
5°^^) for the PSF anisotropy to obtain (C*''°.(5*''°), which should vanish for a perfect PSF 



We show in Figures [T] and [2] results of anisotropic PSF correction for the first HOLICs 
C of spin 1. Figure [1] shows the dimensionless rg\C,\ before (sohd) and after (dashed) the 
anisotropic PSF correction as a function of degree of PSF anisotropy: Aaniso (left panel) and 
6'aniso (right panel). Firstly, it clearly shows that the spin-1 PSF anisotropy induced in mock 
galaxy images is considerably reduced after applying our PSF correction method outhned 
in §31 Secondly, we see a clear trend that the smaller the galaxy size, the more severe the 
anisotropic PSF effect; or that r^lCI increases with decreasing galaxy size, i?gai. In particular, 
when the source size is comparable to the size of PSF (-Rgai = c^iso), the effect is larger about 
one order of magnitude than that of large sources with i?gai = 3(jjso- Figure [2] shows the 
ratio |C*'^°|/|C°'''^| of residual to observed PSF anisotropy in mock galaxy images. We see that, 
overall, the fractional correction factor |C**°|/|C°'''^I is larger for larger galaxy images. Similar 
trends are also found for the second HOLICs 5 of spin 3, as shown in Figures E] and HI 



We apply our flexion analysis method based on the HOLICs moment approach to Subaru 
imaging observations of the cluster A1689. A1689 is a rich cluster of galaxies at a moderately 
low redshift oi z = 0.183, having a large Einstein radius of ~ 45 arcsec {zg ~ 1; Broadhurst 
et al. 2005b). A1689 is one of the best studied lensing clusters (e.g., Tyson & Fisher 1995; 
King, Clowe, & Schneider et al. 2002; Bardeau et al. 2005; Broadhurst et al. 2005a; 




correction. 



5.2. Flexion Analysis of Subaru A1689 Data 
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Broadhurst et al. 2005b; Halkola et al. 2006; Medezinski et al. 2007; Leonard et al. 2007; 
Limousin et al. 2007; Umetsu, Broadhurst, Takada 2007; Umetsu & Broadhurst 2007), and 
therefore serves as an ideal target for testing our flexion analysis pipeline. Deep HST/ACS 
imaging of the central region of A1689 has revealed ~ 100 multiply lensed images of ~ 30 
background galaxies (Broadhurst et al. 2005b), which allowed a detailed reconstruction 
of the mass distribution in the cluster core (10/i~^kpc ^ r ^ 200/i~^kpc). Broadhurst et 
al. (2005a) developed a method for reconstructing the cluster mass profile by combining 
weak-lensing tangential shear and magnification bias measurements, and derived a model- 
independent projected mass profile of the cluster out to its virial radius (r ^ 2h~^ Mpc) based 
on wide-field Subaru/Suprime-Cam data. The combination of weak shear and magnification 
data breaks the mass sheet degeneracy (Broadhurst, Taylor, & Peacock 1995) inherent in 
all reconstruction methods based solely on the shape-distortion information (Schneider & 
Seitz 1995). Broadhurst ct al. (2005a) found that the combined ACS and Subaru profile of 
the cluster is well fitted by an NFW profile (Navarro, Frcnk, & Frenk 1997) with a virial 
mass of Mvir = 1.35 x IO^^H^^Mq and a concentration of c^i^. ~ 13.7, which is significantly 
larger than theoretically expected (cvir — 5) for the standard LCDM model (Bullock ct 
al. 2001; Neto ct al. 2007). Based on these Subaru data, Umetsu & Broadhurst (2007) 
used a maximum likelihood method to reconstruct the two-dimensional mass map of A1689 
from combined shear and magnification data, and found the azimuthally-averaged mass 
profile from the full two-dimensional reconstruction is in good agreement with the earlier 
results from the one- dimensional analysis by Broadhurst et al. (2005a), supporting the 
assumption of quasi- circular symmetry in the projected mass distribution of A1689. Recently, 
Leonard et al. (2007) performed a weak lensing analysis of A1689 based on the ACS data 
by incorporating measurements of fiexion as well as weak shear and strong lensing, and 
their fiexion reconstruction has revealed mass substructures associated with small clumps of 
galaxies, while no fiexion signal has been detected at the cluster center, showing an under 
density around the location of the cD galaxy. 

For our flexion analysis of A1689, we used Suprime-Cam ?''-imaging data (Broadhurst 
et al. 2005a; Umetsu & Broadhurst 2007), covering a field of ~ 30' x 25' with ^.'202 pixel" ^ 
sampling. The seeing FWHM in the co-added i' image is 0'.'88, and the limiting magnitude 
is i' — 25.9 for a Scr detection within a 2" aperture (see Broadhurst et al. 2005a). Since 
the fiexion signal is weaker at larger angular scales (see OUF for detailed discussions), in 
the present fiexion analysis we discarded outer boundaries from the analysis and only used 
the central 3000 x 3000 pixel region, corresponding an angular scale of ~ 10' on a side, or a 
physical scale of 1.3h~^ Mpc at the cluster redshift of 2; = 0.183. 

We used our weak lensing analysis pipehne based on IMCAT (Kaiser et al. 1995) 
extended to include our HOLICs moment method. We selected a stellar sample of — 73 
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objects for measuring the anisotropic/isotropic PSF effects. On the basis of the simulation 
results, we excluded from our background galaxy sample those small objects whose half-light 
radius (r/j) and Gaussian detection radius (r^) are smaller than or comparable to the PSF; 
here we selected galaxies with O'.'G < Vh < 2" and Vg > 0'.'38, whereas the median values 
of stellar and Vg are (vh^) ^ 0'.'48 and (r^*) ^ 0'.'34 using A^* = 73 stars. These lower 
cutoffs in the galaxy size are essential for reliable flexion measurements, because the smaller 
the object, the noisier its shape measurement due to pixelization noise; and the shape of 
an image whose intrinsic size is smaller than or comparable to the PSF size can be highly 
distorted and smeared, as we have seen in §5.11 Faint objects will also yield noisy shape 
measurements, in particular for the case of higher order shape moments. We thus selected 
bright galaxies with 20 < i' < 25 in the AB magnitude system. Further we excluded from 
our background sample those objects whose flexion estimates are significantly larger than 
the model prediction {\F\ ~ 0.1 arcsec"^) using the best-fitting NFW profile derived from 
the joint Subaru and ACS analysis (Broadhurst et al. 2005a); with this model, we set an 
upper cutoff in flexion of |F| < 0.4 arcsec"^, and an upper cutoff in the first HOLICs of 
l^isoj ^ g 02 aj-cggfj-i (^ggg equation [75]). We note that the measurements of spin-3 HOLICs 
6 were found to be quite noisy (see OUF for detailed discussions), so that we discarded the 
second flexion measurements from the present study. Finally, these selection criteria yielded 
a sample of 791 galaxies usable for our flexion analysis, corresponding to a mean surface 
number density of fig = 7.75 arcmin"^. Using the background galaxy sample above we found 
a mean value of (F) = 0.000223 arcsec"^ and a dispersion of ap = 0.11245 arcsec"^. In 
Figure [5] we show the spatial distribution of spin-1 PSF anisotropy as measured from stellar 
images before (left panel) and after (right panel) the anisotropic PSF correction. Figure [6] 
compares the distribution of two components of complex spin-1 PSF anisotropy before (left 
panel) and after(right panel) the anisotropic PSF correction. As clearly shown in Figures O 
and El the higher-order PSF anisotropy in observational data is indeed significant, so that 
one needs to take into account the higher-order PSF anisotropy correction in practical flexion 
measurements. 

We use our first flexion measurements obtained with our moment-based analysis method 
to reconstruct the projected mass distribution of A1689. To do this we utilize the Fourier- 
space relation between the first flexion F and the lensing convergence k (§2.4 of OUF; Bacon 
et al. 2006) with the weak lensing approximation. The field size for the mass reconstruction 
is 9' X 9', sampled with a grid of 256 x 256 pixels, over which the unconstrained k = mode 
is set to zero. Figures [7] and [H] show the E-mode convergence k = ke due to lensing and the 
B-mode convergence kb which is expected to vanish in the weak lensing limit and can thus 
be used to monitor the reconstruction error in the ii^-mode k map. A central 8' x 8' region 
is displayed in Figures [7] and [H The reconstructed k maps were smoothed with a Gaussian 
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filter of FWHM = 0^33. Table [T] lists basic statistics of the reconstructed E- and 5-mode 
K fields measured in the central 8' x 8' field. The rms dispersion in the Gaussian smoothed 
5- mode K map is obtained as as ~ 0.51. The maximum and minimum values in the B- 
mode convergence field are 1.52 and —1.61, corresponding to 2.9cr and —3.1a fluctuations 
(see Table [T]). Figure [7] reveals two signiflcant mass concentrations in the E-mode k map 
associated with clumps of bright galaxies. The first peak has a peak value of ke = 2.66, and 
is detected at 5.2cr significance. This first peak is associated with the central concentration 
of bright cluster galaxies including the cD galaxy, as shown in Figure [9l This central mass 
concentration associated with the brightest cluster galaxies was not detected in the earlier 
ACS fiexion analysis by Leonard et al. (2007). The ACS/Subaru best-fitting NFW model 
predicts k,{zs = 1) ^ 2.5 at ~ O'l. The second peak (4.4o" significance), on the other 
hand, is located ~ 0'.9 to the northeast direction, and is associated with a local clump of 
bright galaxies (see Figure [9]), having a peak value of ke = 2.23. This second mass peak 
has been detected in the earlier lensing studies based on the high resolution HST/ACS data 
(e.g, Broadhurst et al. 2005b; Leonard et al. 2007), and its likely bimodality in the central 
region has been discussed previous studies (Miralda-Escude & Babul 1995; Halkola, Seitz, & 
Pannella 2006; Limousin et al 2007; Saha, Williams, & Ferreras 2007). However, as compared 
to the Subaru weak lensing analysis by Umetsu & Broadhurst (2007), the fiexion-based mass 
reconstruction cannot recover the global cluster structure on larger angular scales ( ^ a few 
arcmin), as demonstrated by OUF. 



6. Discussion and Conclusions 



In the present paper, we have developed a method for weak lensing fiexion analysis by 
fully extending the KSB method to include the measurement of HOLICs (OUF). In particu- 
lar, we take into account explicitly the weight function in calculations of noisy shape moments 
and the effects of spin-1 and spin-3 PSF anisotropics, as well as isotropic PSF smearing, in 
the limit of weak lensing and small PSF anisotropy (g). The higher order weak lensing 
effect induces a centroid shift in the observed image of the background (Goldberg & Bacon 
2005; OUF; Goldberg & Leonard 2007). In weighted moment calculations, this will yield 
in the flexion measurement additional correction terms (relevant to W'{x),W"{x),W"'{x)) 
that must be taken into account by properly expanding the weight function W{x). It is 
found that neglecting these additional terms originated from the Taylor expansion of W{x) 
yields the same result as obtained by Goldberg & Leonard (2007; see Appendix therein). 
We extended the KSB formalism to include the higher-order isotropic and anisotropic PSF 
effects relevant to spin-1 and spin-3 HOLICs by following the prescription given by KSB 
and Bartelmann & Schneider (2001), which provides direct relations between the observable 
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HOLICs and underlying flexion in the weak lensing limit. 

We have implemented in our analysis pipeline our flexion analysis algorithm based on the 
HOLICs moment approach, and tested the reliability and limitation of our PSF correction 
scheme using numerical simulations. Our simulation results show that (i) after applying our 
PSF correction method the PSF-induced anisotropics in HOLICs of mock galaxy images can 
be considerably reduced by a factor of 10-100, depending on the strength of PSF anisotropy, 
(ii) those small galaxies whose angular size is smaller than or comparable to the size of 
PSF suffer from severe anisotropic PSF effects, and that (iii) there is an overall trend that 
the fractional correction factor is larger for larger galaxy images. Therefore, our simulation 
results support the reliability of our PSF-correction scheme and its practical implementation. 

Based on the simulation results, we have applied our flexion analysis pipeline to ground- 
based i' imaging data of the rich cluster A1689 [z = 0.183) taken with Subaru/Suprime- 
Cam. Our flexion analysis of Subaru A1689 data revealed a non-negligible, signiflcant effect 
of higher-order PSF anisotropy induced in stellar images (Figures and E]). It is therefore 
important in practical flexion measurements to quantify and correct for the higher-order 
anisotropic PSF effects. 

Our mass reconstruction from the first-flexion measurements shows two signiflcant (> 
4a) mass structures associated with concentrations of bright galaxies in the central cluster 
region: the flrst peak {5.2a) associated with the central concentration of bright galaxies 
including the cD galaxy, and the second peak (4.4(t) associated with a clump of bright 
galaxies located ~ 1' northeast of the cluster center. This significant detection of the second 
peak confirms earlier ACS results from the strong lensing analysis (Broadhurst et al. 2005b; 
Halkola et al. 2006; Leonard et al. 2007) and the combined strong lensing, weak shear, and 
flexion analysis by Leonard et al. (2007). The central mass peak, however, was not recovered 
in the earlier flexion analysis by Leonard et al. (2007) based on HST/ACS data. Leonard et 
al. (2007) attributed this to their relatively large reconstruction error at the cluster center, 
although they have a very large number density of background galaxies, fig ~ 75 arcmin"^. 
On the other hand, owing to our conservative selection criteria for the background sample, 
the mean number density of background galaxies used for the present analysis is fig = 7.75 
arcmin"^, which is almost one order of magnitude smaller than that of the ACS data, and 
is about 20% — 30% of a typical number density of magnitude/size-selected background 
galaxies usable for the quadrupole shape measurements in ground-based Subaru observations 
{fig ~ 30 — 40arcmin~^). However, we found that it is rather important to remove small/faint 
galaxy images and noisy outliers in flexion measurements since they are likely to be affected 
by the residual PSF anisotropy and/or observational noise in the shape measurement ( §5.ip . 
Besides, the smaller the object, the larger the amplitude of intrinsic flexion contributions. 
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Recall that flexion and HOLICs have a dimension of length inverse: The response to flexion is 
size-dependent, and the amplitude of intrinsic flexion is inversely proportional to the object 
size. Indeed, we find that inclusion of smaller objects results in a noisy reconstruction. 
Similar values of the background number density, Ug ^ lOarcmin"^, have been used in recent 
quadrupole weak lensing analyses based on Subaru observations (e.g., Broadhurst et al. 
2005a; Umetsu & Broadhurst 2007; weak lensing cluster mass measurements of Okabe & 
Umetsu 2008). In their studies only objects redder than the cluster sequence are selected 
in color-magnitude space for their weak lensing analysis, because such a red population is 
expected to comprise only background galaxies {zg ~ 0.9; see, e.g., Medezinski et al. 2007), 
made redder by relatively large /c-corrections and with negligible contamination by cluster 
galaxies (Broadhurst et al. 2005a; Medezinski et al. 2007). However, the smaller number 
of objects implies a coarser angular resolution in the map-making for achieving a proper 
signal-to-noise ratio (e.g., per-pixel S/N ^ 1). With n,, ~ lOarcmin"^ for cluster quadrupole 
weak lensing, typical angular resolutions are about 1 — 2 arcmin (e.g., Gaussian FWHM, 
or boxcar width). On the other hand, flexion measures essentially the gradient of the tidal 
gravitational shear field (i.e., F,G cc (f){r)/r'^), and hence is relatively sensitive to small-scale 
structures. Therefore, our successful reconstruction of the mass substructures with a small 
background density, fig ~ 8arcmin~^, could be attributed to the superior sensitivity of flexion 
to small scale structures (see OUF for detailed discussions) and the here-adopted selection 
criteria for a background galaxy sample for weak lensing flexion analysis. 

Finally, we emphasize that our HOLICs formalism here is different from the earher 
work by Goldberg Sz Leonard (2007) in that (1) additional correction terms for the centroid 
shift, relevant to the derivatives of the weight function, have been included and (2) the 
spin-1 and spin-3 PSF anisotropics, as well as the isotropic PSF smearing, have been taken 
into account under the assumption of small PSF anisotropy (q'[0]). as done in the KSB 
formalism. Our flexion-based mass reconstruction of A1689 demonstrates the power of the 
generalized flexion analysis techniques for quantitative and accurate measurements of the 
weak gravitational lensing effects. 

We thank Masahiro Takada for valuable discussions. We thank the anonymous referee 
for invaluable comments and suggestions. The work is partially supported by the COE 
program at Tohoku University. This work in part supported by the National Science Gouncil 
of Taiwan under the grant NSG95-2112-M-001-074-MY2. 
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A. HOLICs Formalism in the Complex Form 

A.l. Complex Displacement 

In this Appendix, we adopt the following notations to represent complex quantities with 
the degree N and the spin number M: 

= {A9i + iA92)'^{A9,-iA92)'^, (Al) 



= (A/3i + iA/32)^(A/3i + iA/?2) 



N-M 
' 2 



= (A/3i + iA/32)™(A/3i-iA/?2)'^. (A2) 

Here a quantity with the spin number of —M is equivalent to the complex conjugate of 
the corresponding spin-M quantity. Unless otherwise noted, we shall use to represent 
complex quantities in the image plane, and those in the source plane. 

The product of X^ and X^ is expressed as 

vN _ y-N+L 

^M^K — ^M+Ki 

XZX^* = XZt^^ iM>K), 

= (XZ*Xkr^{X^l^^y (M<K). (A3) 

For arbitrary complex numbers W and Z the following identities hold: 

2Re[WX^*]Xk = {WX^* + W*X^)Xk = WX^^^ + W*XZXK, (A4) 
4Re [WXZ*] Re [ZX^*] = [WX^* + W*XZ) {ZX^* + Z*X^) , 

= wz {x^XmT + wz*x^t^^ + w-*zx5+iM + w*z*x^Xm. 

= 2Re[WZ{X^X^^y]+2Re[WZ*X^t^^]. (A5) 



A.2. HOLICs Family 

We summarize in the complex form a family of complex shape moments including 
HOLICs relevant to the weak lensing flexion analysis. 



M = / d^9I{9)W{Xl/a'^) (spin - 0) 
D ^Lmmm^ (spin-l), 
trQ ^ j (f9I{9)W{Xl/a'^)Xl (spin - 0), 



C = i'^'^m^^ f ^l^')^^ (spin - 1) 

^ =/d2^J(^)W^(XoVcT2)Xo4 (spin-0) 

^ ^ (^p.^ _ 2) 

= .f'^'^mw^^l/.^)x^ (3pi^ _ ^) 

vo ^ Jd^0I{0)W{Xl/a^)Xl (spin-0) 

vn = f'i^omwixi/.^)xl ^^p.^ _ 2) 

= (spin - 6) 

rr = /^-^^WH-^(^oVc^-)x7 ^^^^ _ i) 

r,,, = f^^OmWiXi/a^)xr ^^^^ _ 3^ 

ryn ^Jji£immz!m (spin-?). 



A. 3. Differential Operators 

Let us first define the complex gradient operator d as 

Operating d on complex X^ yields: 

/ N+M N-M \ N+M N-M-2 

= d [{Xl)^ {Xl*)^) ^(N-M) {Xl)^ {Xl*)^- = (TV - M)XZl\ 
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N+M 



Af+M — 2 . . ... 

(Ar + M)(X^) 2 {X\*) 2 = (Ar + M)Xi^:l. (A9) 



Similarly, by operating d on the weight function W{XQla^)^ one finds the following: 



9 / 1^2 



(AlO) 



where F and G are the first and the second fiexion defined as F — J^/{1 — k,) and G — 
Q/{1 — k), respectively. Then, the complex displacement in the source plane, A/3, is expressed 
in terms of the lensing convergence, shear and fiexion as 



A/3 = A/3i + iA(32 = ^ {1 - 



Xl - gXl* - ^ {2FX^ + F*Xl + GXl*) 



(All) 



The integration measures in the source and image planes are related in the following way 
d^(5 = (l-fi;)2(^l-2Re[FXii*]-|(7r+^|F|2x2_l|G'|2x2 



+ ^Re [F'^Xl*] -Re [g*FXl] -Re [g*GXl*] —Re [FG*Xl] ) d^O 
fti (1 - nf (1 - 2Re [FXl*] ) d^O, 
to the first order of reduced fiexion (e.g. OUF). 



(A12) 
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Table 1. Basic statistics of the reconstructed E and B mode convergence fields 



field 


mean (/t) 


c 


skewness * 


kurtosis 


minimum 


maximum 


B 


0.058 


0.512 


-0.334 ± 0.2 


0.043 ± 0.4 


-1.61 


1.52 


E 


-0.058 


0.701 


-0.239 ±0.2 


0.125 ±0.4 


-2.09 


2.66 



Note. — The moments are calculated from the convergence within the central 8' x 8' 
region. 

''Skewness defined as ((k — K)^)/a^ . 
''Kurtosis defined as ((k — k)*)/(T* — 3. 
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Fig. 1. — Test of anisotropic PSF correction for the first HOLICs of spin-1, based on 
numerical simulations. No observational noise and lensing signal is included. The values of 
first HOLICs multiplied with the object detection radius, t^ICI; are shown as a function of 
parameters Aaniso (left 6'aniso = 0.92aiso) and 6'aniso {right: Aaniso = 0.1) for a model PSF 
anisotropy. Solid lines indicate the observed values of first HOLICs (C"^^) for Gaussian source 
images smeared with the model PSF, and dashed hnes indicate the residual values (C*''") after 
correcting for the spin-1 PSF anisotropy. No observational noise or lensing signal has been 
added. The PSF consists of an isotropic part P**° described by a truncated Gaussian with 
dispersion aiso and an isotropic part q{9) — Aanisodi/\0\'^ truncated at ^ = Oaniso- Filled 
circles, open triangles, and crosses represent the measurements for a Gaussian source of 
dispersion i?gai = IcTj^o, 2crjso) So-jso respectively. The Cobs is measured with a Gaussian weight 
function of dispersion Vg — R^^t^ + af^^ from the surface brightness distribution smeared 
with the model PSF. 
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Fig. 2. — Ratio of residual to observed spin-1 PSF anisotropy, |C°|/|C°'"'I) ^ function 
of model parameters Aaniso {iefi) and ^^aniso (right) for the PSF anisotropy kernel. The PSF 
consists of an isotropic part P**° described by a truncated Gaussian with dispersion aiso 
and an isotropic part q{9) = Aaniso(^i/\d\'^ truncated at 6' = Oaniso- Filled circles, open 
triangles, and crosses represent the measurements for a Gaussian source of dispersion i?gai = 
laiso, ^aiso, ScTjso, respectively. 
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Fig. 3. — Same as Figure 1 but for the second HOLICs, 6, of spin-3. 
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Fig. 4. — Same as Figure 2 but for the second HOLICs, 6, of spin-3. 
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Fig. 5. — The spin-1 PSF anisotropy field ("^'"^{xj-y) before [left) and after (right) the PSF 
correction over the Subaru z'-band image of A1689. The spin-1 PSF anisotropy was measured 
from stellar shape moments following the HOLlCs formalism outlined in ^ The orientation 
of the vectors shows the direction of the spin-1 anisotropy, and the length is proportional to 
the magnitude of anisotropy. A vector of 0.005pixel~^ is displayed in the inset panel. 
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Fig. 6. — Comparison of spin-1 PSF anisotropy components before {left) and after {right) 
the PSF correction. 
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Fig. 7. — The dimensionless surface mass density k of the galaxy cluster A1689 {z = 
0.183) in the central 4' x 4' region reconstructed using the first flexion observed with Subaru 
telescope/Suprime-Cam. The lowest contour and the contour interval are at a Icr level of the 
reconstruction error 0.51) estimated from the rms of the B-mode reconstruction. The 
black, solid circle in the lower-right corner indicates the Gaussian FWHM (= 0^33) used for 
the mass reconstruction. 
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Fig. 8. — Same as Figure[7]but for the B-mode reconstruction from the first flexion measured 
from Subaru data. The color scale is the same as in Figure [71 The lowest contour and the 
contour interval are at a la level of the B-mode reconstruction (cr^ ^ 0.51) over the 8' x 8' 
region. 
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Fig. 9. — False-color image of the central 4' x 4' cluster region composed of the 
Subaru/Suprime-Cam V and i' images. Overlayed are contours of the lensing K-field re- 
constructed from the first flexion measurements using the z'-band data. The contours are 
spaced in units of 1(t(~ 0.51) reconstruction error estimated from the rms of the i?-mode 
reconstruction. North is to the top, and East to the left. 



